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Abstract 

This paper shows the existence of a periodic orbit with singularity in the sym- 
metric collinear four body problem. In each period of the orbit, there is a binary 
collision (BC) between the inner two bodies and a simultaneous binary collision 
(SBC) of the two clusters on both sides of the origin. The system is regularized 
and the existence is proven by using the implicit function theorem and a continuity 
argument on differential equations of the regularized Hamiltonian. 
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1 Introduction 

In this paper, we prove the existence of a periodic orbit with singularities in the sym- 
metric collinear four body problem of celestial mechanics. The orbit regularly alter- 
nates between two kinds of collisions: a binary collision between the inner two bod- 
ies and a simultaneous binary collision of the two clusters on both sides of the ori- 
gin. Many numerical results related to similar periodic orbits have been published. 
This type of orbit was first discovered in the collinear three-body problem with equal 
mass by Schubart [7| (1956) and subsequently extended to unequal masses by (Henon 
0, 1977) who had also studied the continuation of the family into planar motions 
|1| (Henon, 1976). The unequal mass case was further investigated numerically by 
Mikkola and Hietarinta [4| (1991). There is one such interplay orbit for each mass 
ratio of the collinear three-body problem, which when stable is surrounded by a local 
region of phase space filled with quasiperiodic orbits. These are one of the dominant 
features of the dynamics as illustrated by Hietarinta and Mikkola [3| (1993) . 



Recently, two mathematical proofs of the existence of Schubart orbits have been 
provided: one is a topological proof given by Mockel[5|, 2008 and the other is a vari- 
ational proof given by Venturelli| 10|, 2008. In our paper, a new mathematical proof is 
given for the Schubart-like orbit in the symmetric collinear four body problem. 
Inspired by Sweatman's work[l) and (9), we consider the following symmetric case: 
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the four masses are 1, m, m, and 1 respectively, and they are symmetric about the 
center of the mass, which is assumed to be the origin. By analyzing the regularized 
Hamiltonian, we use the implicit theorem and a continuity argument to show the exis- 
tence of this type of periodic orbit. In section 2.1, we give an important estimate for 
the maximum distance of the outer bodies such that the motion of the inner bodies is 
monotonic. In section 2.2, we introduce the regularized Hamiltonian. In section 3, the 
implicit function theorem is used to show that the net momentum of the cluster at si- 
multaneous binary collision (SBC) is continuous with respect to one particular variable 
R, which is the square root of the maximum distance of the outer bodies. We show that 
there exist distinct values of R such that the net momentum of the cluster at SBC is 
positive for one value of R and negative for the other. Therefore, by continuity there 
exists some R such that the net momentum of the cluster at SBC is 0, which guarantees 
the existence of the periodic orbit because of symmetry. 

2 The Setting and the orbit 

2.1 The setting in Cartesian Coordinate System 

1 m m\ 
— • • • • • * 

—x\ —xi x 2 x\ 

Fig. 1 

From right to left, number the four bodies from 1 to 4. As in the problem setting, 
the masses for the four bodies are 1, m, m,and 1 respectively. The system remains 
symmetrically distributed about the center of mass. The coordinates for the bodies 
are xi, xi, —X2 and — x\ respectively. And body 1 to 4 have velocities x\, X2,—X2,—xi 
respectively. 

The Netowanian equations are 

1 m m 

4x 2 (x\+X2) 2 {x\-X2) 2 

ml 1 

4*2 (xi+x 2 ) 2 (x\-X2) 2 
We will choose the total energy E = — 1 . 

In this paper, we are interested in finding a special periodic orbit with singularities. 
The orbit alternates between binary collision (BC) between the inner two bodies 2 and 
3, and SBC between bodies 1 and 2 and bodies 3 and 4. By introducing a new set of 
transformations, the singularities of BC and SBC in our problem can be regularized. 
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2.2 The Setting in the Regularized System 

We will adopt Sweatman's|9| work to regularize the system.The system has Hamilto- 
nian 

1 o It 1 m 2 2m 2m 



H = — w? H wi 



4 Am 2x\ 2x2 X1+X2 x\ — X2 ' 

where w\ — 2x\ and W2 — 2m±2 are the conjugate momenta to x\ and X2- Introduce a 
canonical transformation 

q l =x l -x 2 , q 2 =2x 2 , pi=w u p 2 = ^(wi + w 2 ). 

This results in a new form for the Hamiltonian 

1 , p? P\P? Po 2m m 2 2m 1 
H= l + -)!-L- 1 lF2 i F2 



tn 4 m m q\ q2 qi+q2 2q\+q2 

To regularize the equations of motion, we introduce a Levi-Civita type of canonical 
transformation 

Qi = ^ql, P,=2Q, Pi (i = 1,2), 

and we also replace time t by the new independent variable s which is given by j- — 
q\q2- In the extended phase space, this produces a regularized Hamiltonian 

dt , 

_ 1 2ja , QlPl-*QiQ2PiPi+4Q\Pl 2 2 2 2mQ 2 Q 2 Q 2 Q 2 2 2 
16 QlPl + 16m mQl 2mQ2 ~ Q 2 + Q 2 ~ 2Q\ + Q 2 ~ QlQlE > 

where E = — 1 is the total energy. 

We start at BC with initial conditions 

xi (0) = A, x 2 (0)=0, *i(0)=0, i 2 (0) = +oo 
which is a singular point. 

To analyze the motion, it is necessary to deal with the singularity in the regularized 
coordinate system. The corresponding initial conditions at s = in this new coordinate 
system are: 

Qi(0)=R, G 2 (0) = 0, Pi(0) = 0, ft(0) = 2m3, 
where A = R 2 . 



3 



2.3 Estimation of A 



Intuitively, if A is big enough, there will be multiple BCs before the first SBC happens. 
In order to find the desired orbit, we will have to give an estimation of A such that there 
will be no BC for t £ (0,fi], where t\ is the time of first SBC. 

Definition: Let a body has velocity v = at time t* . If there exists an interval [t a ,tb], 
such that t a < t* < tj,, and v is positive for t E [?«,?*) and is negative for t S (t* ,tb], or 
v is negative for t E [t a ,t*) and is positive for t E (t* ,tb], then we call t* is the turning 
time and the position of the body at t* is called the turning point. 



-X2 



Turning Point 



A'2 



Theorem 1. There exists an Aq, such that the second body has no turning point for 
t E (0,?i] whenever < A < Aq, where t\ is the time when the first SBC happens. 
Further, the second body will have at least one turning point for t E (0,?i] if A > Aq. 



Proof. Consider the Newtonian equations of x\ and X2'. 



x\ 



X2 = 



1 



Ax\ {x\+X2) 2 {x\-X2) 2 



4x1 (x\+X2) 2 (x\— X2) 2 



(1) 



(2) 



In order to get an upper bound of A, we consider a necessary for the second body 
changing direction. Assume t = t* < t\ is the time when ±2 = 0, then x% < for 
tE (0,f*]. 

Letxi(f*) = ax2(t*) with a > 1. Because X2(t*) < 0, a must satisfy 

1 1 



m 
4 



i.e. 



(a + 1) 2 (a -I) 2 
\6a<m(a 2 -\) 2 . 



<0 



(3) 

Also,forf E [0,t*], xl/x2>xi(t*)/x2{t*) =a. By the setting of the problem, a > 1. 
Rewrite equation[T]as: 



-xi 



4x? {X1+X2) 2 (X1-X2) 2 



1 2m(l+|) 



(1-|)2 



4 



Note 2m(l + y) /( 1 — y) 2 is increasing with respect to y if y < 1 . In our case, y - 
x\jx\ < 1/a 2 < 1, then 



2m(l 



r 2 

4) 



2^2 



< 



1 , 2(! + ^) 



4 ' (l-^) 2 



Therefore, 



-i'l < -=■ 



1 2ma 2 (a 2 + l) 
4 + (a 2 -l) 2 

1 2a 2 (a 2 +l) 



(a 2 -!) 2 



Let jti(O) = A and x\(t*) =A\ < A. Since xi is decreasing for t € [0,f*], i.e. < 
as r e [0,f*]. The following inequality is true for t € [0,?*]: 



1 

x\ 



1 2wa 2 (a 2 + l) 
4 + (a 2 -l) 2 



> x\x\, 



integrate this from t = to t = t * to get: 

1 2ma 2 (a 2 + l) 
4 + (a 2 -l) 2 



' / 1 \ f 

— x\— j at > x\x\dt, 



Then 



As E = -1, then 



(1-1) 
% A 



(1-1) 
y Ai A' 



1 2ma 2 (a 2 + l) 
4 + (a 2 -l) 2 

1 4wa 2 (a 2 + l) 



> fun- 



(a 2 -l) 2 



-i=if(0- 
=*?(**)- 



1 w 2m 2m 

2xi 2x2 x\+X2 x\—xi 



1 

A? 



1 

a7 



By inequality |U 



K(---) 
~ V Ai A ; 



1 ( 4ma 2 (a 2 + l) 
2 



1 

a7 



(fl 2 -l) 2 



(fl 2 


-I) 2 


m 2 a 


4ma 




a 2 - 



2m 


2m 

i 


a 


i 


_ l 


4ma 2 

i 




a 2 


-i 




1 




2 

m- 


"aT 


2 + 


~2 



fl 2 -l 

1 4ma 2 (a 2 + l) 
2 + (a 2 -l) 2 



(4) 
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1< 



Then 



1 ma[16a — m(a 2 — l) 2 ] 1 
= T x 2(a 2 -l) 2 A 

Applying inequality [3] to this gives, 
1 ma [16a — m(a 2 — l) 2 ] 1 



1 4ma 2 (a 2 +l) 

2 + (a 2 -!) 2 



2(a 2 -l) 2 



1 4ma 2 (a 2 + l) 

2 + (a 2 -!) 2 



< 



1 4ma 2 (a 2 + l) 
2 + (a 2 -l) 2 



1 4ma 2 (a 2 + l) 
A "2 + (a 2 -l) 2 ■ 



(5) 



By inequalities |3] and |5] 



1 12m 

A > - + 4m - 



8m 



'--1 {a 2 -!) 2 ' 

where a 4 — 2a 2 — + 1 > 0. In other words, if the following two inequalities A < 
l/2+4m+12m/(a 2 - 1) + 8m/(a 2 - l) 2 and a 4 - 2a 2 - \6a/m + 1 < hold, x 2 will 
be monotonically increasing as? £ (0,fi]. Choose Aq, such that X2{t*) =X2(t*) =0for 
some time f* 6 (0,fi], that is 



a —2a — ■ 



16a 



1=0 



m 



(6) 



and 



8m 



1 12m 
Ao>-+4m + - I - T - 1 (fl2 _ i)2 . 



When A — Aq, derivative the equation |2] with respect to t evaluated at t*: 

„. = mx 2 (t*) 2(x x (t*)+x 2 {t*)) _ 2(x l (t*)~x 2 (t*)) 
X 24(f*) [xi(t*)+x 2 (t*)} 3 [x^n-x^t*)} 3 ' 

Since x 2 {t*)=0, x 2 (t*) = 0, xi(t*) < 0, xi(f*) > x 2 (f*) > 0, 

1 1 



"i2(0 = 2*i(0 



[x!(f*)+x 2 (f*)] 3 [xi(f*)-x 2 (f*)] 3 



>0. 



Hence, when A = Ao, jc 2 will increase monotonically for f £ (0, f i ] . Therefore, x 2 is 
increasing for? £ (0,fi] wheneverA < Ao, where t\ is the time of the first SBC, andx 2 
will have at least one local maximum for t £ (0, fo) when A > Ao- □ 

Remark: In the special case m = 1, equation [6] becomes a 4 — 2a 2 — 16a + 1 = 
with a > 1, a « 2.766. Then 

1 4a 2 (a 2 + l) 
A >- + T2 - TF i> 6.485. 
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3 Existence of the Periodic Orbit 



Recall that in Section l2l2l 

where E = — 1 is the total energy. The initial conditions at s = are: 

Ql(P)=R, g 2 (0)=0, P\ (0)=0, P 2 (0)=2mi 

By Theorem [T] when < /? = < \/Ao, 22 = V2~X2 is increasing from s = to 
s = s\, where s\ is the time when the first SBC happens. 

The equations of motion from the regularized Hamiltonian Y are: 

Q[ = l -^QlPi - ^-QiQiPi, (7) 

8 m Am 

Q2 = ^-Q 2 iP2-^-QiQiP2, (8) 

2m 4m 

.1 l AmQ\(A 2Q 1 Q$ , 

^ ^ 5^ - ^+ 2 ^ + T^fgi + 7^%! - 22.6!, (9) 

where ' is the derivative with respect to s. The initial conditions are 

Gi(o)=fl, (22(0) =o, Pi(o)=o, p 2 (o)=2. 

At the time Si when the first SBC happens, 

&(*) = <), C22(ji)=J?i, P l (s 1 )=-^=. 

\/2m + 2 

To prove the existence of the periodic orbit, we are going to find a value of R, such 
that P 2 (si) =0 when gi(si) =0. 

Theorem 2. For f/ze regularized Hamiltonian T, let the initial condition be P\ (0) = 0, 
Pz(0) =2, 2i(0) = /? > anof 22(0) = 0, where R e (0, V^o]- Assume at time s\ 
\ i (/?), Q2 > 0, Qi = 0, which means SBC occurs. Also, assume Q2 > Ofor < 5 < Si . 
77zen P2(s\^R) is a continuous function ofR. 
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Proof. Since the Hamiltonian Y is regularized, the solution P, = Pj(s,R) and Q\ = 
Qi(s,R) are continuous functions with respect to s and R. We are going to show 
s\ = s\ (R) is a continuous function of R. In order to apply the implicit function theo- 
rem for Qi = Qi(si,R) = 0, we need to show that (dQi/ds)(si,R) ^ 0. 
By the regularized Hamiltonian T, 



-fa- \( S uR)- l Pi \(-H,R) 



1 + m 



1 



m 4m 



si M) 



Note that for fixed R, T = at any time s. At s = si, Q\ = Qi(si,R) = 0, then 
Pi =Pi(s u R) =4. Therefore, 



"37 ^ 



1 + ra , 1 



m+ 1 



Q 2 2 (s u R)>0. 



By the implicit function theorem, si is a continuous function of R. Then P2(si,R) 
is also a continuous function of R. □ 

Corollary 3. There exists R such that ftCsi) = Pi{s\ ,R) = 0. 

Proof. First, we are going to show that there exists an P > such that ) > 0. 

From equations ITlfTOl 

(PiQi +P2Qi)' = P[Qi +PiQ\ +P2Q2+P2Q2 



-4mQ 2 2 + 2m 2 Q 2 + 2Q 2 Q 2 2 



2m 



1 



-2 



Note that A > jci (f ) >x 2 (t) > for? g [0,fi], thenO < gi < P = VA and 2*i(fi) = 
2x 2 (fi) = Q\{s\ )=R\<2A= 2R 2 , that is < Q 2 < Ri < V2R for s e [0,si]. Thus 



2m 



> 



2m 



Q\ + Q 2 2 2Q 2 + Q 2 2 ~3R 2 4R 2 



Choose R — 



2m 



Q\ + Ql 2Q 2 + Q 2 



4m 



so (P1Q1 +P2Q2)' > for s G [0, Ji], or P1Q1 +P2Q2 is increasing for s € [0,si]. 
From the initial conditions, 



(P1Q1+P2Q2) U=o= 0. 



Hence, 



o< (Pie 1 +p 2 e 2 ) i« 1 =*ift(jo 

Therefore, when R = , P 2 (*i) > . 
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Next, we will show P 2 (si) is negative when R 2 =Aq. 



At A = Aq, by the proof of theoremQ] there exists a time t* <t\, such that ±2 (t * ) = 
andii(f*) < 0. Thenii (f*) +x 2 (t*) < 0. Consider the sum of the Newtonian equations 
□ andEl 

1 1 2 

X\ +X2 = T T — -, rr < 0, 

4x 2 ( Xl+X2 )2 

which means x\(t) + x 2 {i) is a decreasing function with respect to t. Hence, + 
x%{h) <xx(t*)+x 2 {t*)<Q. Notefhatft(*i)/[2(22(ji)]=p2(ri)=ii(»i)+jfe(»i)<0, 
and g 2 (*l) >0, soP 2 (*l) < 0. 

Therefore, by continuity, there must exist an R, such that ^2(^1 ) =0 where s\ is the 
time when the first SBC happens. 

□ 

Theorem 4. If R satisfies P2(si,R) — 0, then the orbit will be a one -dimensional 
Schubart-like periodic orbit. 

Proof. At time s = 0, a BC happens between bodies 2 and 3. At time s — s\, a SBC 
occurs. Since the system is regularized, the solution {P{,Qi} (i = 1,2) is continuous. 



At time s = 0, 

Qi(0)=R, (22(0) =0, Pi(0)=0, P 2 (0)=2mi 
At time s=si, 

C2i(ii)=0, Qi(si)=Ri, P l (si) = --^L=, P 2 (si) = 0, 

\/2m + 2 

where R\ is a positive number. From the Hamiltonian T, we can see that Q[(s\) — 
\Qi{QiP\ - Gift) < 0, Q' 2 (si) = \Qi(2QiP 2 - Q 2 P\) = 0. This means that Q 2 {s x ) 
is the local maximum of Q.2- In other words, when time passes s\, Q2 will decrease. 
Similarly, Q\ will decrease when time passes s\. 
At the time s = 52 when the second SBC occurs, 

Qi{s 2 ) = -R, Qi(s2)=0, Pi{s 2 )=0, P 2 (s 2 ) = -2,J. 

Compare the motion for s £ [0,si] and the motion for s £ [51,52]- By the uniqueness of 
the regularized Hamiltonian system and symmetry, the orbit from 5 = 51 to 5 = 52 will 
be the same trajectory from 5 = 51 to 5 = by reversing the direction of each velocity. 
Further, 52 = 25 1. By symmetry and uniqueness again, at time 5 = 35 1, 

£2i(3ji) = 0, Q2(3ji) = -*i, ft(3*i) = ft(3*i) = 0. 

V 2m + 2 

At time 5 = 4s\, 

Q l {4s 1 )=R, Q 2 {4si)=0, Pt(4ji) = 0, P 2 {4s l )=2,J, 
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which is exactly the same as the initial condition at s = 0. Then the orbit from s = to 
s = 4si will generate one period. 

□ 

The following figure is a picture of the periodic solution for m — 1 in terms of 
{Qi,Q2,P\,Pi}- The initial conditions are 

fii(0) =2.295, &(()) = 0, Pi(0) = 0, ft(0)=2. 



■ Ql 

02 



3 - i \ I 




-5 1 1 1 1 1 1 

2 4 6 8 10 12 

Figure 1 : The horizontal axis represents time s. 
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